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ABSTRACT 


Because of the desire to improve operational readiness 
and to Simultaneously reduce support costs, there 1S a 
great deal of interest in the military Services in imple- 


menting multi-echelon models for determination of adequate 


but economical stocking levels for spare parts. Two models 
presently used--METRIC and MOD-METRICQ--are inefficient and 
require excessive time for computation. In an attempt to 


solve these deficiencies a heuristic model was developed at 
the Naval Postgraduate School. The main purpose of this 
thesis is to compare the characteristics and performance of 
the simple heuristic model with the MTERIC solution. 
Comparisons with METRIC revealed that the heuristic 
model was much more efficitnt computationally, but the solu- 
tion was frequently far inferior to that obtained by METRIC. 
The comparison indicate strongly that base(shipboard) stock 
levels as determined by existing allowance models are larger 


than are needed in an integrated system. 
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I. INTRODUCTION 


A. BACKGROUND 


"No matter how large our forces or how modern our mili- 
tary equipment, if our forces are not ready to fight, or if 
they cannot be sustained once engaged, we have no real 
combat capability" [Ref. 1]. These words were used in the 
Secretary of Defense's annual report to Congress(FY83) to 
promote the idea that a balance must be maintained force 
modernization and support of existing forces [Ref. 2]. The 
Support of existing forces is an isSue of high interest to 
military managers. Continually they try to find the most 
efficient way to allocate the budget that is assigned for 
Spare parts. Many inventory models have been built to help 
determine the proper stockage levels for recoverable items 
which maximize performance subject to a given spares invest- 
ment. One of. these models is METRIC(Multi-Echelon Technique 
for Recoverable Items Control). 

METRIC is a mathematical model translated into a 
computer program, capable of determining base and depot 
stock levels for a #group of recoverable items. Its 
governing purpose is to optimize system performance for 
specified levels of system investment. METRIC is designed 
for application at the weapon system level, where a partic- 
ular line item may be demanded at several bases and the 
bases are supported by one central depot [Ref. 3: p. 123]. 

Another model for spare parts allocation is MOD-METRIC 
which is a model for a multi-item, two-echelon two-indenture 
inventory system. MOD=eMETRIG as an extension of METRIC, 
which permits the explicit consideration of a hierarchical 


Berts structure. The hierarchical parts structure means 


that the recoverable items themselves contain recoverable 
items. . 

The objectives of the MOD-METRIC model are to describe 
the logistic relationship between the components and _ the 
final assembly, and to compute base and depot spare stock 
levels for all items with explicit consideration of this 
logistics relationship [Ret jeep. 472). The difference 
between the METRIC and MOD-METRIC models is the manner in 
which the average resupply time is expressed. Chapter II of 
this paper provides a detailed review,of these two models. 

Both models are excellent tools for finding efficient 
Spare parts stockage solutions for small-sized problems. 
However, because they uSe a recursive solution technique, 
each run takes much computation time, making the models 
infeasible to use for large-sized problems. To overcome 
these computational difficulties, and to provide a workable 
solution for practical problems an heuristic stockage model 
has been developed at the Naval Postgraduate School 
[Ref. 5]. 


Be JPUREOSE 


The main purpose of this paper is to compare the solu- 
tions given by the simple heuristic N.P.S. model to the 
optimal solutions given by the METRIC. The comparison of 
the models will look at both efficiency in terms of computa- 


tional time and an the quality of the solutions obtained. 


C. PREVIEW 


Chapter II of this paper provides a detailed review of 
the METRIC and MOD-METRIC models. Included in that chapter 
are the assumptions of the models and the solution tech- 


niques used for each. 


LO 


Chapter III addresses the problem of making the METRIC 
and the N.P.S. models more efficient by approximating the 
expressions for base backorder days and by estimating base 
stock levels using simple multiple regression equation 
instead of the recurSive computations required for the exact 
solutions. The regression equations uSe as inputs only the 
ready rate, and the demand rate. 

Chapter IV describes the assumptions and solution proce- 
dure of the heuristic N.P.S. model for solving the same 
stockage allocation problem. E 

Chapter V presents the results of example computer runs 
ich illustrate how close the approximate solutions of the 
heuristic N.P.S. model are to the optimal METRIC solution. 
The comparisons include a modification of the N.P.S. solu- 
muom Usine marginal analysis to incorporate unit costs. 


Chapter VI summarizes the thesis and presents. 


‘lee 


IIT. DISCUSSION OF METRIC AND MOD=MniRae 
This chapter describes the METRIC and MOD-METRIC inven- 
tory models. For each model the following subjects are 
discussed: maintenance system structure, assumptions, data 
requirements and notation, the objective function and solu- 
tion ‘techanaaqies ME ER EG ees more fully described by 
Sherbrooke [Ref. 3], and MOD-METRIC by Muckstadt [Ref. 4]. 


A. METRIC 
1. Maintenance System Structure 


Consider “the multi-echelon maintenance structure as 
used in the METRIC model. In a multi-echelon structure 
stocking/maintenance facilities are organized in a hierarch- 
ical structure according to supply/maintenance flows which 
are represented as an arborescent network. 

When a unit fails at base level there is a prob- 
ability r that it can be repaired at the base, and a prob- 
ability l-r that it must be returned to the depot. Because 
of the typical high costs and low demand for items, the 
inventory stockage policy is (S-1,S), which means that items 
are not batched for repair or resupply request. If there is 
inventory available at the base, a serviceable item replaces 
the failed item. If no inventory is available, the equip- 
ment will be inoperable until anitem is repaired. in 
elther case, the failed item immediately begins base repair; 
or 1S sent to the depot. When the failed item is sent to 
the depot a resupply request to the depot is issued. If the 
depot has a serviceable part, it will send the item to the 
base immediately. If the depot does not have the item 


avallable, it will send a serviceable unit to the base after 


IZ. 


a repair is accomplished. ies sprocessais shown in figure 
Zz, 1 


Base Base 
---) component Component 
Failed Repair Inventory 
Component 
Remove & 


Repair 


4 
oom em @ fe ® @ ow 


Depot Depot 
Component Component 
Repair Inventory 


Faitled Component 
Serviceable Component 





Figure 2.1] METRIC Repair Process. 


The depot and base stock level of aniitem will 
affect the performance of the system. As stock levels 
increase the average resupply time of the failed item will 
decrease. The objective of the METRIC inodel is to determine 
the base and depot stock levels of every item for a given 
budget constraint such that the total backorder delay at the 


bases is minimized. 
2. Data Requirements and Notations 


The METRIC model requires several input data items 
for implementation. The required data are listed below; 
m = The number of recoverable items. 


n = Number of bases. 


©?) 
i 


ihe costvot, 1tem i. 


r= Probability that a failure of item i at base j can 


be repaired by base; it follows that l - Ti is 


ae, 


Note; 


probability that item i cannot be repaired at base j. 
Expected order and ship time for item i from 
base j to the depot(days). 


Expected base repair time of item i at base j(days). 
Expected depot repair time of item i(days). 


Expected number of demands for item i at location 
j(demand/day) j = 0)1(,2 23. 

Stock level for item 1 at location i 

j =O nl eee | _ 


Expected resupply time for item i at base j. 
Backorders for item i at location j ; j] = O}10232200 


subscript i is used to index the items and subscript 
j refers to the different bases. The subscript j=0 
refers to the depot. 


Assumptions 


A stationary compound Poisson probability distribu- 
tion describes the demand process for each item(for 
our comparisons later we consider only the case in 
which demands are standard Poisson distributed). 
There is no lateral resupply between bases. 

There are no condemnations(all failed parts are 
repaired). 

A failure of one type of item is statistically inde- 
pendent of those that occur for any other type of 
1tem. | 

Repair times are statistically independent. 


There iS no waiting or batching of items’ before 


repair is started on an item(infinite channel queuing 


assumption). 
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7. The level at which repair is performed depends only 
on the complexity of the repair. 

8. Items and bases may have different military essen- 
tialities; however, items are normally considered to 
be equally essential. [Ref. 4: p. 474] 


4. Objective Function 


The objective function used by METRIC is to minimize 
the sum of expected backorder days on all recoverable items 
at all bases pertinent to a specific weapon system. . A base 
Pie-omcer Lor am | item exists any Cime there is an unsatis- 
fied demand for that item at the base level. Depot back- 
orders are of interest only insofar as they affect base 
backorders. 

To derive’ the ‘expression for expected backorder days 
the following mathematical background is needed. 

A key result for both METRIC and MOD-METRIC concerns 
mien PLObability distribution of the number of units in 
resupply. For the case where demands are Poisson distrib- 
uted, Palm's Theorem gives this distribution under’ the 
assumptions stated for the METRIC and MOD-METRIC model. 

Palm's theorem(stated in the context of the stockage 


problem); 


Let s be the spare stock for an item where demands are 
Poisson distributed with customer arrival rate i. Let 
ib be the probability distribution of resupply time 
rE: h mean T. Then, in the backorder case, with and 
S-1,S)  stockage olicy Ehe steady-state probability of 
fumes if resupply 1s Poisson distributed with param- 
eter AT; l.e. 


AT} exp(-AT 
We.) )-spix units in resupply) = Cat} exp(-AT) 


SS OO 


Patcemeeneo wm aisttibutaon of the number of units in 


resupply depends on the resupply time only through its mean 
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value. The actual resupply time distribution as) moe 
required. 

Feeney and Sherbrooke [Ref. 6] extended Palm's 
Theorem for the case where demands are compound Poisson 
distributed under the assumption that all demands placed by 
a given customer have the same resupply time. They show, in 
this case, that the resulting distribution of the number of 
units in resupply iS compound Poisson with parameter AT. 

in addition. Feeney and Sherbrooke looked at the 
special case in which the demands per,customer are logarith- 
mically distributed with variancee-to-mean (7 aememar They 
enoned that for this special case the resulting distribution 
for the number of units in resupply is negative binomial 


With parameters q and k = AT/I1n q; i.e. 


p(x units in resupply) = 


k+x-1)! (q-1l 
(ere-1)t (ae go... g > 


Using the result given by Palm's Theorem, we can 
then compute easily the expression for the steady state 
expected number of backorders. Let S be the number of units 
of stock allocated to a base and let T be the mean resupply 
time for the base. The number of backorders at a given time 
will then be zero if the demand during the resupply time is 


less than or equal to S and the number of backorders will be 


(x-S) if the demand is larger than S. The expected number 
of backorders at any given time, denoted B(S,T), is then 
given by; 
oe) 
B(S,T) = > (x-S)p(x|AT) (eqn 2.1) 
x=st+l 


Since the value of My for a base is assumed known, 
we need to compute only the mean resupply time Tij for item 


1 at base 7: 
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feces nk tl - rs N*(A + expected depot delay time). 
1) iL 3 1d ud 


Except for the expected depot delay time all the 
variables are assumed known. . The depot delay time is zero 
if the depot has infinite stock. If the depot has no stock, 
the time is D(depot repair time). If the depot has finite 
stock S then the expected depot backorders is 


ee) 
= = ae eqn 2.2 
Peg 41% > Sig) PCI A; Ps) alias: 
i0 
: n 
where dA. = GL se, ee 
oO, j= 1 ah A) 


Equation 2.2 can be interpreted as depot backorder 
days per day [Ref: 7]. Thus , when we divide this number by 
depot demand per day(rA. ) we get the expected depot back- 


order days per demand: 


oO 
i Me eae). oe 
She: Oy 


_ om tL 
10 10 
Sherbrooke argues that this average delay will be some 
fraction of the depot repair time D and uses the notation 
Bios, )D for this expression to emphasize this fact. And 


so the resupply time can be expressed as follows; 


T (S.)=r R_ + (1-r,_)(A.. + 8,(8,_)D.). (eqn 2.4) 
a) 6|CULO sy Meelis | 1j 1j a OS 


Now we can write the equation for the expected 
number of backorders for item i at base j when the depot has 


Baecond the base has $|; stock levels; 


N] 


B..(S_.;S 


See) Pe (x| x TT... (S..)). (Sane!) 
mam ij iLO K=S4j sk i 1 wee 


A 


As mentioned above, the objective of METRIC ismee 
minimize the sum of backorders for all items i and for all 
bases j within a budget constraint. The METRIC problems is 
igi Peemmems acto, S is greater than or equal to zeroufor alia 


ebalel 9) | \sieulic.o. 


n 
Minimize > yal B . (Saee So eee (eqn 2.6) 
T=] 7J=1 15 1) 3xe 


_ m n 
Subject to CAs < B : 
: ran = te la} 


Where B is the available budget. 


5. Solution Technique 


The METRIC problem can be solved by using either 
marginal allocation or the generalized Lagrangian Multiplier 
method [Ref. 3: p. 1S Srl We describe the generalized 
Lagrange multiplier method since that is what is suggested 
by the authors of METRIC. Rirsts we state the theorem by 
Everett which is the basis for the generalized Lagrange 
multiplier method. 

Everett theorem [Ref. 8]: 


Let S bea aoe eginip Leto rs SO BEET | of possible 
strategies or actions and let H(x) be the ay off VGae 
Ded ey whieh acerues cere ee the stra egy X Ea 

5 E » a 


Let C(x) be the resource required by s BECES 8 nd let 
C be the maximum amount of the resource available. We want 
to 

Maximize H(x): for all x€ §S (eqn 2.7) 


Subject to Co). Ge 


This problem can be expressed as an _unconstrained problem 
for a given Lagrangian multiplier 6> 0 as follows: 


Maximize H(x) - @C(x). (eqn 2.8) 


Ibis 


If an x can be found which maximizes equation 2.8, then this 
5 aS also the solution for the constrained problem, equation 


Everett s Eheorem is useful for solving the METRIC 
problem. According to Everett's Theorem we can express our 
original problem, equation 2.6 as follows; 

Main = ae 
i=l j=l 
where @6@< Q.0 


Sa (eqn 2.9) 
Tag 


es Go ‘3 Oo = 
Sere, ed 


c 
sa aa =) 


a 


Different choices of @ lead to different resource levels, 
and it is necessary to adjust them by trial and error to 
achieve a given constraint. Therefore we need to’ solve 
problem 2.9 for several values of the multiplier @. Neale 
value which provides the resource level that ius closest to 


the budget constraint will be chosen. 


Since our problem is separable in the items, problem 
2.9 can be solved for each item separately. That is, we can 


solve the m subproblems, 


Ni 
oven (6) 4S5..) - Cs) -Se.s_). (eqn 2.10) 
jek ia ele 1.0 1, a Wg 1 10 


The solution technique is outlined below: 


Step 1. In the description g refers to the trial value 
for the Lagrange multiplier at step k. 


Select a starting value 8 for the multiplier; 


where @ < 0. A reasonable choice for a will 
be presented later. 

peep 2. Establish an upper bound on S_.. This will be 
presented later also. ; 

Step 3. Given S__ and 9 , determine the base stock level 
Belg Seek Desa (SJ. may start at zero). Since 
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B. (S..,8._) is diseretely convex for vaseivenm. . 
Tj Ly 


(see [Ref. 9: p. 260] ): 


B(S...8,..) -@C S&S <= B (S) +15.) epee eee 
ak Ej 


ig a © jena ag] i)... 1 10 ie 
or 
B. (8S, .+1,8,.) = By Gee ec eee (eqn 2.11) 
ane e el) 10 19) Se k i 


To find the optimum stock level for each base, 
given depot stock Sjg and a increase a from 
O to the smallest integer which satisfies 
inequality 2.11 


Step 4. Increase depot stock oa by one unit, and return 
to step 3. Continue increasing depot stock until 
it reaches the upper bound established in step 2. 

Step 5. Choose optimal S__ for each item. The optimum 
S UJ aiidetts oor: scoonuie base stock levels are 
those values which minimize problem 2.12 


nN 


Min B oS - @CS§ Smee) On es) eqn 2.12 
& | iy! ij io? kag Kea io? ied ) 


Step 6. Change the item and go to Step 2. 


Step /. Compute the required total cost for buying all 
S for a given @. If total cost 159022 ceen as 
ea (where «is an acceptable prespeemiredad aes 
ference total cost and budget), then stop. The 
current Lagrange multiplier gives optimal 
solution for those resources actually required by 
the solution SS... Otherwise choose a new @_ and 
SO —EOSteD o.. eee con search procedure should 


be used to determine a new value for 


BeeGnorce of a Multiplier and Depot Stock Level Range 


Sie sauniGtonmeat equation 2.5 into equation 2.11 


gives 


> (x - (8..+ 1) p(xfaT(s,_)} 
Lj +0 


x=S 
1J 
€O 
ee eC pix (Ar(S )} 2 6 C, 
X=5S_ ag: i0 k 1 
1J 
Which reduce to : 
oe) 
eS ee 8, 
a P(xJAT(S 9) 2 06; 
wey) 
Add 1.0 to each side gives 
oO .. 
(Ome AT(S )} 2 1.0 +04.c. 
pee) kd 
1J 
Finally, we rewrite the left side: 
“ 
mas LCS yp ee 0s tee Orme eqne2. 1s 
page N80? ki - 


If we define the "ready rate" to be probability that 
the quantity of an item demanded during a resupply time is 
less than or equal to the stock level Sq; , we see from equa- 
tion 2.13 that the METRIC solution forces the ready rate for 


each item to be at least as large as 1.0 + 6C,. Therefore, 
meewe choose a4 minimum ready rate, a lower bound on 6 can 
be estimated from inequality 2.13 . This minimum value of 


can be used as the starting Lagrange multiplier in step 1 of 
the solution technique. 

We also need to establish a range of values for 
depot stock level Sjg . Because the depot demand rate is 
and repair time is Dj , the average number of units of items 


mieaepot resupply is A,D;. Empirical evidence in the runs 
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which are discussed in Chapter V suggests  thabms seme oe 
never be larger than the greatest integer less jGiane ee 
accomplish mean supply response time goals in the neighbor- 
Hood NOme melas. 


The minimum should be taken to be QO, and the search 
conducted over the values O07192,.-4, (ia to find Bie 
METRIC “sot aom- 


B. MOD-METRIC 
1. Maintenance System Structure 


Consider a maintenance system consisting of two 
echelons--depot and bases. Assume that an end item is 
divided into several repair modules. For example, if an end 
item iS an aircraft engine it may have, modules for intake, 
combustion and exhaust. If an engine fails at a base, it is 
replaced by a serviceable engine from base stock. The 
failed engine then goes immediatly to either base repair or 
is shipped to the depot and a resupply request is sent from 
the baseto the depot. 

When the engine iS repaired at the base, it is 
assumed that one of the modules is faulty. A serviceable 
module from the base stock, if available, will replace the 
failed module, and the repaired engine is placed in base 
engine stock. The failed module that is removed from the 
engine at a base is repaired at the base or is sent to the 
depot anda resupply request is submitted to the depot. 
This entire repair process is shown in figure 2.2 

Now, there are several different stock levels which 
affect system performance. Since engine and modules have 
different functions in the repair process, each stock level 


affects the system performance differently. 


ae 


Base Engine Re- 
air, Remove 
eplace Module 


Legend 


Failed Engine 
Repaired Engine 
Failed Module 
—-—-—-— Repaired Module 





Rrgure) 2a MOD-METRIC Repair Process. 


2. Data Requirements and Notation 


The data items required by MOD-METRIC are listed 


below. 


m = Number of modules associate with the end item. The 
subscript i will index-modules, i = 0 represents the 
end item; 

n = Number of bases. The subscript j will index bases, 


j = O represents the depot; 
A .=Average number of daily removals of module i at base j; 
1J 


r = Probability a failure of module i at base j requires 
1Jj 
base repair; 


ZS 


R. = Average base repair time for module i at base j; 


1j 
A .= Average order and ship time for module i at base j 
ee EO Ene depot, 
EEG Average depot repair time of module i; 
a. = Cost of module i MG. means cost of engine); 
ae Stock level for module i at location j,; j = OJ, ee 
eg Average resupply time for module i at base j. 


a 


5. , ASSUMpeEdons 


The assumptions stated in METRIC are also applicable 
to MOD-METRIC except andy ee Instead of compound Poisson 
demand in METRIC, MOD-METRIC assumes that the demand process 
is the simple Poisson process. METRIC assumes that each 
item has the same military essentiality. In MOD-METRIC (Gia 
essentialities of end items and modules are explicitly 
expressed through the equation that represents each item's 
contribution to the end item's resupply time. Furthermore, 
MOD-METRIC assumes that if an engine requires repair and 
that repair is made at the base level, the probability that 


more than one module requires repair iS zero. 


4. Obiective Function 


AS in METRIC, the objective is to minimize the tora! 


expected engine's (or end item) base backorder days summed 
over all bases. As before, the minimization is subject towa 
budget constraint. In deriving the objective function only 


one end item and its modules are considered. 
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If we can find the expected base resupply time of 
the end item at base j, then we can use equation 2.1 for 
the base backorders of the end item. An expression for T 


is.as follows; 


aT r R_ + (A_. + expected depot delay time). 
Oj Os 07 0j 


When the depot has Spo , the expected depot delay 
time per demand is very similar to equation 2.3 . Let 
co) D = Expected depot delay per demand. 
00° 00 | ; 
inven, 
> 
8(S D = --- SS x|A__D : eqn 2.14 
(Soy Pog NT Ses 00 Pil ngPog?* ~— Sed 


O00 O00 


- 


nN 
where = > (l-r )A is the expected depot demand. 
00 a= Oy. 207 


Average base repair time for end item at base j, Ro}; 
is equal to the average remove and replace time, given the 
necessary module is available, plus the expected delay due 


to the unavailability of the module which is required to 


repair the engine. Therefore Ro; = Boj + Ao; where 
ar = the average repair time at base i if modules are 
J 
available; 
eae = the average delay in base engine repair due to 
J 


the unavailability of a needed module. 
Let the expected delay in engine base repair time due to a 
backorder on module i at base j be represented by a:- 


\J 
Then 


OO 
i. - 2 Cmts) p(x A T..) (eqn 2.15) 
ale} eae te 1 1j jest deg 


ED) 


T= feR (2) Cee (eqn 2.16) 
4 aes fenaiey ij ij 

i 

ee x “Sao A, .D 

07 Se a1 0 eo 


and 


Th 
r = iD pe ws r 
10 a! iy ij 


a 


The-expected delay in engine repair at base j due to modules 


s 
™m 
ee ee a ee (éqn 2.17) 
10 YO es) 
Oj Oj 


Thus we have shown that the average resupply time for an 


engine, Toj » can be expressed as: 


ik = r B + A tp 2qait Pe Ales 
0j o;$ 0j 03? ( q ) 


Now, the problem of MOD-METRIC is to find Si; > Q shich can 


be expressed mathematically. 


TN oc 
Minimize Ea pa (x - So JPG) (eqme2s 17m) 
joa sos Oj oie) 
J 


mn Tn 
Subject to an = CS. Sap 
L= 
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39- oolution Teehnique 


Pew (Soliton, technique suggested by Muckstadt 
[Ref. 4] is outlined in this section. 


Problem 2.19 is mot separable because Toj is a 


complex function of the Sij : The approach taken by 
Muckstadt was to partition the problem og ETe: two 
subproblems--the module subproblem and the end item 


subproblem. The solution algorithm is as follows; 


Step 1. Set up minimum investment levels for modules and 
end items. Let those be g and £, respectively. 
Let z be the total expected backorders for the 
end item when the base stock levels for the end 
item are S_.. Set a budget increment for modules, 
be tse Re fgg trial value for total system modules 


iivestmeame and let zZ =o . 


Step 2. Solve the module subproblem given the minimum 
budget g = g. That is, find Sij 2 0 which 
Minimize > 3 S  (x-S__) P(xla,.T..(S__)) 
inimize x~ Scie 
‘gh to ae i) eee a0 
18), 


- 
Subject to Gs < ae 
= j=0 i ag 


Step 3. Solve the end item subproblem. Find oe QO as 
ineede MEBRTEG solution which 


n Ce) 
Minimize z => 2 Ca cence As Rees.) ) 
Ie Fe 0 j 0j Oj ° 00 
a 


nh 
Subject to Es 205 eee 
fo DOH : 


et 


where T_ is calculated using the module stock levels 
ee | 

determined in step 2 and given S__. Compute also the 

value of z, the value of z associated with this opti- 


mal “Se aWwe lon. 


Step Geet ece Zz, eo directly to steppe 
Otherwise let z = z and retain the corresponding 
stock levels as the incumbent stock levels. Then 
go to step 5, 

Step 5. Increase ¢ by) b Deeb a’ < f, then terminate. 
Otherwise return to step 2. If the algorithm seam 
on this step, the optimal stock levels and the 
associated minimum expected backorders will be 
those saved from step 4. In step 2 an Gptimizatuom 
problem in solved in which a portion of the budget 
& , is allocated among the modules to determine 
depot and base module stock levels. The objective 
of this subproblem is to determine that division 
of e’ which minimize the sum of the expected delay 
function. The optimal valuegotme. is found = by 
trial and error by searching through the 


integers. 


Note that the form of the suboptimization problem in 
step 2 is exactly the same as the METRIC problem. Thus , 
step 2 can be solved using the techniques that were used in 
ME TRIG. 

The engine delay time due to unavailable modules at 
a base(4o; ) is fixed by solving step 2 for a trial modus 
budget. This means that the expected engine repair time at 
a base dis fixed “in problem 2. lowe. Thus step 3 is also 
exactly same as the METRIC problem for only one item. 
Consequently, the METRIC solution technique can be used in 


Step 3 again. 
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TII. REGRESSION TECHNIQUE 


As observed in the METRIC and MOD-METRIC solution tech- 
niques, a lot of calculations are required to compute B (Sij 
moig ) and to find Sij (Gain nie), 0,1 ,2...n). Because 
of the computation time it is difficult to determine the 
Beares allocation for thousands of items. It would be much 
faster if the computationally heavy ,recursive calculations 
feeoe replaced by a mom-recursive procedure. ‘indes 
chapter presents resuits of an effort to approximate the 


exact computations uSing multiple regression equations. 


For estimation of the regression equations the following 
procedure was used( the i,j subscripts are suppressed since 


the results apply to each base and item separately): 


Ime oelect p= AT(S). 
Vary the base stock level over the integers; S=0 1,2,.. 
3. For each value of S, find the true ready rate, 


RR(s) = (xsS|#), and backorder function 


S (x-S)p(x|#). 
x 


=.) 
4. Select those values of S for which the ready rate is in 
a predetermined interval like, say 80% to 95%. 


5. Repeat steps 1-4 a range of values of uz. 


The set of values selected at step 4 constituted the regre- 
ssion data base for the predetermined ready rate subinter- 
mae. NOw, with this data base, MINITAB was used for two 


sevarate regression analyses. 


mie sase stock, S, was regressed on »# and ready rate RR. 
2). Base backorders, B, was regressed on i ere Sa sie 
AMG . 


Ze 


Tables I and II show the regression output from MINITAB. 
They show how each variable contributes to the prediction of 
S and eBe 


TABLE I 
Regression Equation for Base Stock 
(5. D0<h<7 (5 Sago oR Rae 
THE REGRESSION “EQUATIONS: 
Ye=2- 11.3 + 21:23 3) =e 
Y iS a expected base stock 


COEF 


ante 
ge 


NT 
oy: 
CH 
7 





Separate equations were estimated for arbitary selected 
Subintervals for wand RR. We selected 5 subintervals forpu 
and 5 subintervals for RR giving a total of 25 possible sets 


of equations. 


Table III shows how accurately the regression equations 
estimate the actual § and B (S,#,RR ) for a given u , and 
10 

Table III reads as follows; if a base has 7 spare parts 
and its » is 5.15, then its ready rate is 854 and Dackordgen 
days) dageme 0529 77 ~ A regression equation is also used to 
estimate the base stock and backorder days assuming an 85% 


ready rate and the »w from the first column of Table III 


‘Table III uses the equation that appeared at table I 
and table III . 
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TABLE LI 
Regression Equation for Backorders 


ooo Rio) 


(oe Q<H<7.5, 
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Table III shows that the regression equation gives very good 


estimates of the values for both the base stock and the 


backorder days. Thus, if we are able to use such equations 
in solving multi echelon inventory problems, we can reduce 
the computation time very much. The N.P.S. model does use 


these equations and which will be presented at the next 


chapter. 
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iveenEURISTICG MODEL 


A. BACKGROUND 


As we have seen in the chapter II and III, the METRIC 
and MOD-METRIC solutions require extensive recursive Como 
tation. For example, the first recursive process involves a 
search for the optimal base stock level given the depot 
stock level and a value for the Lagrangian multiplier. This 
process is repeated every time the depot stock is changed--a 
Second recursion. Furthermore, the above two recursive 
processes are repeated for each item and each value of the 
Lagrangian multiplier.. Because of these recursions, much 
computational time is required, and it could be difficult to 
find the optimal solution for a system having thousands of 
items(the typical aircraft has about 2,000 work unit coded 
repairable item [Ref. 10] ). 

Recall that in the METRIC model it was shown that the 
@peimal base stock level was the smallest integer value of Sj 
which satisfies 


Pex Ss 5 hs 1 + HC .. 
1J 1J a. 


If we ignore the integrality of demand we can interpret this 
inequality as providing the same ready rate for each base. 
Now, this has the altinative interpretation that given a 
value of the Lagrange multiplier and a fixed depot stock 
level, the METRIC solution corresponds to the determination 
of the largest ready rate which is budget feasible. 

Now, aside from minimizing base backorder Gays, itis 
frequently the case that a base may wish to achieve a speci- 
fied ready rate assuming no Support from a higher 
echelon(depot). Such is the case, for example, with ship- 
board(base) allowance list models which stock to provide a 


90% ready rate for a 90 day period of time. 


S)6: 


The N.P.S. stockage model incorporates this notion vores 
Specified base ready rate and attempts to determine the 
depot stock level Sg required to achieve amean supply 
response time(MSRT) goal. For our study, the goal wae 
selected arbitrarily to be 125 hours (5.2 daycur 


B. SOLUTION TECHNIQUE 


The base ready rate is a function of the amount of depot 
stock(the depot stock determines the Tij ) and the amount of 
base stock. Also the optimal base stock is a functionwen 
the-ready rate and depot. stock. Furthermore, as shown in 
equation 2.4, the base backorder is a function of Sij and 
TAs CS jo ire Therefore, the first step of the heuristic model 
is to determine the base stock as a function of fixed depam 
Stock and ready rate. The second step is to determine the 
value of the base backorder days as a function of the fixed 
depot stock and the base stock which was found in the first 
step. 

The necessary values can be found recursively using the 
Same expressions used in METRIC or they can be approximated 
very accurately and rapidly using regression functions. The 
N.P.S. model saves much computation time when solving for 
the optimal base stock levels and base backorders by using 
the regression equations described in the previous Chapter. 

The depot stock leves for a given item by beginning at 0 
and incrementing depot stock by one unit until the average 
base MSRT first reaches the specified MSRT goal. Thus S,,18 


the smallest integer such that 


| 
Be < MSRT 
sr Ba5 8057800? 4 2a; 2 


Where B_ (S ane is th total expected base backorder days 
1jJ, OJ | 
for item i at base j when base j ha S units and the depot 


1J 
has ae units of stock. MSRT is the specifiedsiskl s2oaam 
a & 
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This procedure is repeated for each item one-by-one. 

The big difference between the METRIC and N.P.S. models 
is that METRIC minimizing the expected number of base back- 
order days for a fixed budget, whereas the N.P.S. model 
attempts to find the minimum depot stock levels required to 
Mewtevewa Gesired supply response time goal for each item. 
A comparison of the METRIC and N.P.S. models is presented in 
Smapter 5. i 


oe 6 ALGORITHM . 


a 


-Below is an outline of the solution technique used by 
the N.P.S. model. PEiec eNO PreCrerence 15 given tO any item 
by the N.P.S. model and since there is no budget constraint, 
no consideration need be given to the unit cost. © Therefore 


the same solution procedure is repeated for each item. 


Step 1. Find the smallest non-negative integers S_, 
1J 
to). 2 eee, Stich that 

Ries os...) 2 059 
al 
where T = 90 days. the regression equations 


described in chapter III can be used for this 


Step. Set the depot stock S equal to QO. 


Step 2. Determine pte and BME oe) 16) 9) ee Le ees 0 We 
(The regression equations can be used to approxi- 
mate the backorder functions required by this 


step). 


Step 3. Compute the demand-weighted average MSRT across 


the n bases. 
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Steeped 


If MSRT is -less than 125 hours then Step ene cums 
rent S_ is the desired solution. If MSRT greater 
than 125 hours,” amepeace Si by one unit and go to 
step 2. 


Repeat for each item. 
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V. COMPARISON OF THE MODELS 


inweemdereta evalliatre the heuristic, non-optimal.,...N.P.S. 
model we selected several sample data sets and computed the 
N.P.S. solution for each. We then repeated the computations 
using the METRIC solution procedure. Since the METRIC solu- 
tion is optimal, we can evaluate the quality of the N.P.S. 
selution by comparing its performance, to that of the MEIREE 
mode | . The results of several comparisons are contained in 
this chapter. 

Example problems consist of three items which are to be 
stocked at one depot and three bases. For the purpose of 
mielustration, four data sets are constructed. First the 
N.P.S. model solves the problem and computes the associate 
budget, then this budget is used aS a constraint in the 
METRIC model. 


A. DATA SET 1: LONG ORDER-SHIP TIME AND LOW PROBABILITY OF 
BASE REPAIR 


Table IV gives the complete input data for the first 
data set. The first example considers long order-ship times 
and zero probability of base repair. Thus, all repairs must 
take place at the depot. Item 1 has a low demand rate, item 
2 has © a medium demand rate and item 3 has a high demand 
rate. The demand rates are consistent through the four 
sample data sets. 

Recall from chapter II the following notations: 


r = probability that item i can be repaired at base j, 
1J 

meee ene expected repair time in days for item i at base j, 
1Jj 

A. .= the average order and ship time in days from base j 


1j 
to the depot for item i, 


o7 


TABLE IV 
INPUT DATA 1 


[tem Base 


ie 


NIEBDS QOH ANE 
COfD NWeH NOL 
QO@Q@ @Q@@ 7 Oe@ ie 


0.0 
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= 


D = the average depot repair time in days for item i, 


i 
A = the failure rate(failures per day) of item i at 
base j. 


Tables V and VI give the N.P.S. and METRI@G@resu tee 
Observe that there are significant differences in the 
Stockage allocations between the two models. However, the 
difference in performance(MSRT) is not so large. Notice 
also the large difference in computational times; One 
seconds for N.P.S. solution and over 20 times as long(5.8/7 
seconds) for the METRIC solution. 


B. DATA SET 2: LONG ORDER-SHIP TIME AND HIGH PROBABILI TT: aan 
BASE REPAIR 


Table VII gives the data for the second test data set. 
Characteristics of this data set are long ordering and ship- 
ping times and a high probability of base repair. Thus, the 
bases will rarely need to rely on the depot for Support. 
Intuitively, very little stock will be needed at the depot 


level. 


a0 


TABLE V 
Results of N.P.S. Model for Data 1 


system Average MSRT 
Time for Computer Run 


TABLE VI 
Pe cowlinswGi MeaewleumMecder. for Data lL 


+ 
ihetal Cost): 
system Average MSRT 
Time for Computer Run 





As expected, both solutions place little stock at the 
depot. The N.P.S. model puts zero stock there for each 
item. Observe also that the system average MSRT values are 
much lower than the response time goal of 5.2 days(125 
Mmaurs } . This is because of the self-supporting capability 
of the bases. As before, the N.P.S. solution required a 


small fraction of the time required by the METRIC solution. 
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- TABLES ee 
INPUT DATA 2 


H 
a 


Item Base 


OOo @O2OGQ  @eeaqg 
POPE Hr OOO 
SIRI QOWRH AUS 
COP NOR NOL 
OQo®d Oaoqa eae 
CONJCO \OCOCO ~wr.\oce 
OD OOM OOM 
WILD POPs Pohopye 
Unio COo~po Wn 
OoOood OOO DOO 
Ooo QoQ @Weqg@ 
OoOooO DOO DOO 


9 
9 
9 
9 
9 
9 
9 
9 
9 


TABLES V ELE 
Results of N.P.S. Medel for epacae 


system Average MSR 
Time for Computer Run 





C. DATA SET 3: SHORT ORDER-SHIP TIME AND ZEROGPROG Cilla 
OF BASE REPAIR 


Table X gives the data for the third sample data set. 
This set is characterized by short ordering and shipping 
times and low probabilities of base repair. Table XI and 
XII are the N.P.S. and METRIC model results. 
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TABLE 1X 
Recmltcmonm MErkre Model for Data 2 


Item 2 
Item 3 
+ 
Total Cost : : 
System Average MSRT : ; days 
Time for Computer Run : : econds 


a 


TABLE X 
INPUT DATA 3 


Item base 


Pee 
Ky 
a 
> 


GOe CQOeG7 eee 
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As with the previous data set we see that the N.P.S. 


model gives zero stock at the depot. METRIC, on the other 
hand, carries less stock at the bases and positive stock 
levels at the depot. Both models easily satify the MSRT 


meat of 125 hours. 
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TABLE XL 
Results of N.P.S. Medel tor Data 


System Average MSRT 
Time for Computer Run 


TABLES one 
Results of METRIC Model for Data 3 


+ 

Total Cost : . O- 
System Average MSRT : ; days 
Time for Computer Run : : econds 





The results from analyses of the last two data sets show 
that the system mean supply response times can be made 
Significantly smaller than the 125 hour goal, even with zero 
Stock at the depot. This is because the stockage levels at 
the bases are very high. This example points out a short- 
coming of the N.P.S. model(and perhaps of the current Navy 


allowance list models); too much stock is positioned at the 
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Saaseemeenbe  jostatwed by an MSRP coal of 125 hours.* As a 
result of the pe eee vac ons), the N.P.S. model was modi- 
fied for the last data test case by reducing the base ready 
fate from 904 to 7/024. 


ieee DATA SET 4: ZERO BASE REPAIR AND LONG DEPOT REPAIR TIME 


Table XIII gives the data for the fourth test data set. 
Characteristics of this data set are short ordering and 
shipping times, zero probability of base repair and long 


a 


depot repair time. 


TAD Ee kt 
INPUT DATA 4 


H 
a 
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A comparison of the results above to those obtained for 
Giese set 3 shows that the MSRT goal of 125 hours(5.2 days) 
can be achieved even with longer depot repair times’ and 
lower base ready rates(/0% versus 90%) at a reduced invest- 


ment cost. This supports the argument above that the base 


*This comment is contingent on the validit of the 
ene ytons made by METRIC model. Dpepartiveu lar, the 
(S=155)) ordering policy assumption and the Poisson demand 
assumption. 
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TABLE wy 
Results of N.P.S. Model for Baca 


System Average MSR: 
Time for Computer Run 


TABU Eee 
Results of METRIC Models tor pata 4 


+ 
EOE AL COS.tag, 
oystem Average MSRT : 
Time for Computer Run 





stock levels are too high. Even with the reduction in base 
ready rate to 70%, the solutions above show that the METRIC 
solution puts less stock at the bases and more at the depot 
than doés the WN.P.S. sotutror. The solution obtained by 
METRIC is nearly 20% better than the N.P.S. solutions 
Thus, additional modification of the N.P.S. model to reduge 


base stock even more should be considered. 
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Suservemineene MEtRrG Solution that the base stock allo- 
eapwons ton, item by are identical to those inthe N.P.S. 
Sorlut ion. Since the METRIC Lagrange multiplier solution 
gives each base the same ready rate for a given item it is 
clear that the ready rate selected by METRIC for item 1 is 
70%, the same as that used by the N.P.S. model for each 
item. However, METRIC selected ready rates for items 2 and 
3 which were smaller than 70%(the base stockage levels for 
items 2 and 3 are smaller than those given by the N.P.S. 
Sepeuit lon) . These smaller ready rates are a reflection of 
the higher costs for items 2 and 3. Thus, the METRIC solu- 
tion illustrates that the ready rates selected for the base 
stock levels in the N.P.S. model should not necessarily be 
the same for each item, but should be a function of the unit 
costs. The cheapér items should get higher ready rates than 


the more expensive items. 


Pee ODIFIED N.P.S. MODEL 


moe areindlwAnalytic Solution 


The sample analyses discussed in the previous 
Section suggested that the N.P.S. heuristic model could 
perhaps be improved if the item mean supply times were 
allowed to vary depending on unit costs. Thus, it might be 
better to provide greater protection for the less expensive 
items and less protection for the more expensive items. The 
METRIC solution does discriminate in this Pash yon. 
attempting to provide the greatest performance per dollar 
invested. Therefore, in an effort to improve the efficiency 
of the N.P.S. model and to incorporate consideration of the 
unit costs, the N.P.S. model was modified by using marginal 
analysis to determine the depot stockage levels. 

The modified model determined the base stock levels 


just as before. Then, to determine the depot stock levels 
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for multiple items an iterative process was used to select 
sequentially that item at each step which provides’ the 
greatest reduction in total base backorder days per dollar 
invested. The necessary dollars are allocated to purchase 
one unit of the selected item and the process is repeated 


until the overall MSRT objective was achieved. 


2. Mathematical Description 


Let (d, ,dy »-++5,dy)be the depot stock levels for 
items 1,2, caymeeemespectivelys (Note: This is a change in 
the notation used previously to emphasize that the only 
decision variables incorporated in the marginal analytic 
solution are the depot stock levels for the n items). Let 
B, (dj) be the total base backorder days for item i whengd 
units are stocked at the depot. Recall that the depot stock 
level influences the base backorder expression only through 
the depot resupply time Ti: To emphasize this, we express 
the depot resupply time for item i and base j as Thi (dive 
The total expected base backorder days for item i can then 


be Wrltten Aas: 


N 
B(d.)= 2 > (x-S__)P. (xla,.T..(d.)). (ean 5-1) 
1 Vea j=1-x=5 0 ay Tae pt fea: Wig Je ana 


1j 


Now, we have seen previously that the average depot 
delay per demand in satisfying base j when the depot has dj 


Units of stock 1s piven by. 


T (4) =A. + (3 Ged) BGI Dee 
1j 1 a) x=d_ a aly at a: 
BE 


sh 
where A. = 24 (l-r ) A... is the total expected demand at 
i j= ij ij 
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the depot per unit time for item i; D; is the average depot 
repair time for item i; and Ajj is the order and ship time 


from the depot to base j. 


Define 4B (d) = B (d_) - B (d +l) 
ak iG 1 alk al al, 


to be the reduction in total base backorder delay for item 1 
Miethe depot stock level is increased from d units to d,+l 
units. Because the impact of a change in depot stock on the 
total base backorders is felt only indirectly through a 
reduction in Tij (di), no useful analytical simplification of 
Bey Cd; ) could be determined. Nevertheless, it can be 
computed easily directly. 
Finally, let A= = > 


§2 
demand over all. items and bases. We are mow ready to 


Mj be the total expected 


describe the marginal analytic procedure. 


Step 1. Determine the base stock levels S._. as before 
1J 
from ready rate considerations. Set the initial 
depot stockage vector to be (0,0,...,0); i.e. 


Gi emo meodoea = 1,25... .m. 
at 
maepe 2. secompute T (d +1) for i=1,2,..,m aac) a= le 2 ea 
iy. 1 


eeep 3. Compute B (d )/C for i=1,2,..,m and let k be 


a etaasa ll at 
PitmetmaGdex tot whieh this ratio iS maximum. 


Step 4. Let a = a, ale. 


m 
Step 5. Compute MSRT = (om Bh el i 
ele “Seis oe 


Step 6. If MSRT < Goal, stop. Otherwise go to step 2. 


3. Example of the Marginal Analysis Procedure 


The marginal analysis modification was applied to 
the data set 1 examined in the previous section. The 


resulting allocation are presented in Table XVI 
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TABLE XVI 
Results of the Marginal Analysis 


System Average MSRT : 
Time for Computer Run 
“NP. oO. Mmogels =: RI =.32 
“METRIC MODEL :°RSRi 3-24. 





As may seen . from this Table, the marginal analysis 
modification provided only a small improvement over the 
N.P.S. model(MSRT decreased from 5.011 days to 4.939 days 
and costs were reduced by $900) and performance still falls 
Significantly short of the MEIRUTG sotursen, Observe also, 
as expected, that the marginal analysis soltuutton aga 
increase the stockage levels of the less expensive items and 
decreased the stockage level of the most expensive item. 
This 1S evident from the item MSRT values. While ite 
clear that any solution for which budgets and unit costswaa. 
considered should improve overall effectiveness in a budget 
constrained environment, it may not be desirable to allow 
the more expensive items to suffer in terms of stockage 
Support. Incorporation of a workable essentiality coding 
policy could be used to override the impact observed above. 
While the increased performance obtained by using 
the marginal analysis modification was not significant, keep 
in mind that the marginal allocation was applied only to the 
depot stock levels. It is reasonable to assume that addi- 
tional improvement could be obtained with similar modifica- 


tion to the base stock allocation procedure. 
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VI. CONCLUSIONS 


There were basically three major objectives of this 
thesis. The first was to review the literature and _ to 
understand the existing models for controlling the stockage 
decisions in multi-item, multi-echelon inventory system for 
repairable items. The second objective was to develop the 
computer programs needed to implement , the most promising of 
Mdemcxisting algorithms. The last objective was to evaluate 
fhe heuristic N.P.S. model by comparing it to the baseline 
solution given by the existing algorithm. 

The review of the literature revealed that basically all 
of the existing models are derived from the METRIC family 
developed at RAND in 1968. Therefore, the METRIC model was 
selected as our baseline. That algorithm was programmed and 
is operational on the N.P.S. IBM 370 model 3033 computer. 
Appendix A of this thesis contains a listing of the PL/I 
source code for that program. The N.P.S. model proposed by 
Apple [Ref. 5] was also programmed-in PL/I on the N.P.S. 
computer system and is included in this thesis as Appendix 
B.. 

Execution of the METRIC and N.P.S. models revealed 
extensive recursion in the computations of stockage levels 
and the backorder expression. Since it is well known that 
one of the major problems preventing widespread acceptance 
of the METRIC family of models is the extensive computa- 
tional time required to solve for the optimal solution for 
reasonably-sized weapon systems, we sought to determine if 
accurate approximations could be used in place of much of 
the recursion. 

Multiple regression equations were developed for estima- 


tion of the expected base backorder days as a function of 
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the mean demand, the stockage level, and the ready rau 
Regression equation were also developed for estimation of 
the base stock levels, asa function of the mean demand and 
the ready rate. The accuracy of both sets of equations was 
Very Shae ny This suggests that significant reductions 
CONpuULer et lMewe in elie METRIC model(or others)could be 
achieved by replacing the exact recursive computations with 
the approximate regression equations. 

In order to accomplish the third object ee several 


Sample data sets were created for use in testing the METRIC 


and IN. Fas." mecels, Four of the data sets are described in 
Ghapter V. Analysis of the results of the sample data sets 
revealed some shortcomings of the proposed N.P.S. model. 


One of those shortcomings, complete disregard for item costs 
was corrected by modifying the N.P.S. model to include costs 
in a marginal analysis procedure for determining depot stock 
allocations. This modification achieved only minor improve- 
ment but Since it was restricted to the depot stock deci- 
Sions it requires additionalwseudy. Analysis also revealed 
that the base stock levels are consistently too high in the 
N.P.S. model. Significant reductions in the base ready rate 
from 90% to 70% and less were made with the support system 
Still able to achieve specified mean supply response time 
goals at significantly reduced costs. We also observed that 
the ready rate protection period, nominally stated as 90 
days seems longer than necessary with the types of values 
used for order and ship times and repair times in our exam- 
ples. Clearly, both the ready rate and the length of the 
protection period should be a function of the other input 
parameters of the model, like order and ship times, base 
Fepair WeLEobabal Tey. and base and depot repair times. 
Selection of a uniform "90% protection for 90 days “pole 
for determining shipboard stock will sometimes over protect 


or under protect significantly. This might have dramatic 
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consequences on presently used allowance list models used by 
the military services. 

An evaluation of the computation times required for 
determining the N.P.S. and METRIC solutions showed the times 
for the METRIC solution to be as much as 20 times as great 
as those required for the N.P.S. Souliitaren,. Thus, the 
N.P.S. solution does hold promise for implementation in the 
real world with large weapon systems. 

In general, the performance of the N.P.S. model was 
somewhat disappointing in comparison to the METRIC solution. 
The results of this thesis indicate that the N.P.S. model 


requires additional work in the areas mentioned above. 


OL 


APPENDIX A 
COMPUTER PROGRAM FOR METRIC MODEL BY PL/I 


*PpL/C ATR SOURCE L=9000 P=150 T=(1,30)XREF 
BBO: PROCEDURE OPTIONS (MAIN); 
DCL(COST(15,10,1),BREP(10,10,2) | BDEM( 10) 1033) Bee c eee 
bord(10,10,4) ,bPROB(10,10,5) ,DREP( 10710) Goes 
DCL(DDEM,UB(20),RS,LS,QK,P(70) , SPACE(20) , BRST(20) , DDELY, 
SUMCOST , BODD,RR(30,6,45,3),MSRT)FLOAT; 
DCL(PP,PPP,M,WIJ(30, 6,45, 3),U,KO,LAM(70),PX,POISSON, 
BUDGET , QK1,SUM5 , BOD1, SUM6 , READY _RATE(15) ) FLOAT; 
DCL(Y3,SUM1.,SUM2 ,MU(20) ,BOD2 ,QK2,TOTAL_ DEMAND ) FLOAT 
ne gia ie al ab @a)410))) 
DCL(1I,J,K,L,%,Y,Z2,V,F,HAN, LAMBDA, NUMBER) EEE 
DCL (Si0.S1) , 1D, BASE) FIXED itce 
DCL(SUMCOST1,SUMCOST2 (40) , SPACE1,SUM3 ,SUM4,BOD,W1,W(40) ) 
Elost; 
DCL (MMSRT,TOTAL AVG MSRT,TOTAL MSRT,SUM DEMAND) FLOAT; 


LAMBDA = 11; 

OK -=<(0,6-1.0)7 200. 
ED} = 352 

BASE =3;3 

GET: Dist A(BUDGED): 


DOT = ee Or LD; 
DOG 1916 -BASE, 
je a 
GET LIST(COST(1,J,K),BREP(1,J,K+)) .eDEM eee 
BORD(I,J,K+3), BPROB(1,J,K+4 eS DREP Ci eee 
ENP: 
END; 
SUMCOST -=s00:: 


a2 


DO WHILE (ABS(BUDGET -SUMCOST)>749.0); 
LAM (LAMBDA) =QK; 


PUleSKkrP EDIT (OR) (X(1).F(11L¢8); 
SUMeEOSTL =0.0. 
SUMGOsST =0.0; 
TOTAL DEMAND =0.0; 
howe MoRT = 0.0; 
SUM DEMAND =0.0; 
DOSHAN = -1 TO 2D; . a 
TOTAL DEMAND = 0.0; 
DDEM =0; 
DOG =- Ler Owe Aor: 
PP =BDEM(HAN,J,3) * (1-BPROB(HAN,J,5)); 
DDEM = DDEM +PP; 
TOTAL DEMAND = TOTAL DEMAND + BDEM(HAN,J,3); 
END; 
U =DDEM *DREP(HAN,1,6); 


SPACE1] = 1000.0; 


POR Weiner GhO «< 18) 
PX = EXP(-U); 
SUMZ a= Os oO: 
12 SILOM 
THEN DO; 
Dee = 0210 SLO= 1. 
SUM ese (SiO - X) = PX: 
rae 0 = Pees (X+1); 
END; 
END; 


DDELY =SUM2 /DDEM; 


SUR ee= 0 70': 
SUM4=0.0; 
SUM5 = 0.0; 


Deo 


RS = QK * COST(HAN,1,1); 

BOD=0.0; 

DO I = 1 TO BASE; 
SPACE(I) =0.0; 

END; 


DO G = 1 TO ;BaoE: 
ES 10) 
DO Y =l TO BASE; 
M =BREP(HAN,Y,2) *“BPROB(HAN,Y,5) + (1.0- 
bprob(han,y,5))*(BORD(HAN,Y,4)+DDELY); 
BRST(Y) = M; 
END; 
PUT SKIP(1) LIST('THE BASE MU'); 
_DQ Z = .1 TO BASE; 
UB(Z) = BDEM(HAN,Z,3) *BRST(Z); 
END; 


Si =.0 : 
he lie 
¥3uq1 205 


POISSON = EXP(-UB(G)); 
PROB SUM = POISSON; 
RR(LAMBDA,HAN,SI0+1,G) = PROB SUM; 
DO WHILE (PROB SUM < 1.0+RS); 

Y3 = Y3 *UB(G) / K; 

PROB = POISSON *Y3; 

PROB SUM = PROB SUM + PROB; 

Ke ke 

SIJ = SIJ +l; 
END; 

RR(LAMBDA,HAN,SI0O+1,G) = PROB SUM; 


WIJ(LAMBDA,HAN,SI0+1,G) = SJ; 
PX = EXP ([UBCGne 
SPACE(G) = UB(G) Seman 
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is SILI 218, 
THEN] DO ; 
DOPX=0 10 SII =1; 
SEAGIG )) -SPACE(G) + (SiJ e2eX )*PX; 
PX = UB(G) * PX/(X+1); 
END. 
END: 


SUM5 =SUM5 + SPACE(G); 

BOD = BOD+SPACE(G) -QK*COST(HAN,1,1)*SIJ; 
END: | ; 
DO Z=1 TO BASE; 

PU Ge splen (UB (Z mec yk (9.5) )i; 
END; 
EO enous) Cex Gl); ECE, ) ); 
MMSRT = SUMS / TOTAL DEMAND; 
G =1; 
BODD=BOD - QK*COST(HAN,1,1)*SIO; 
PUT EDIT (MMSRT,BODD)(X(2),F(9,5),X(2),F(9,5)); 
IF BODD<SPACE1 
THEN DO; 
DO I = 1 TO BASE; 
READY RATE (1) = 
RR(LAMBDA ,HAN,SI0+1,1); 
END; 
BOD1 = SUMS; 
SPACE] =BODD; 
POs T© BASE: 
BBOD(1I) =SPACE(I); 
END; 
Wie=Si0: 
NUMBER =0.0; 
DOmmE-mEeTO BASE = 
W(1I)=WIJ (LAMBDA ,HAN,SI0+1,1); 
NUMBER=W(1)+NUMBER; 


3) 


EN: 

END; 
oLO -—-S10y 
END; 


SUMCOST1=COST(HAN,1,1)*(NUMBER+W1); 
BODZ = BOD2Z + BODI; 
PUT SKIP EDIT(HAN,WL)(X(10) , F(Z), X(2 eee 
DO I = 1 TC BASE. 

PUT EDIT(W(I),READY RATE(IL)) 

(X(2) F(Z) eae 

END 
MSRT = BOD1/TOTAL DEMAND; 
PUT EDIT (MSRIZSoUMeceS TS 

(X(2) FG, 55k eer 
SUMCOST =SUMCOST = —sUMees i. 
SLO. =0 205 
SUM_DEMAND = SUM_DEMAND + TOTAL DEMAND; 
TOTAL MSRT TOTAL Moki sion we TOTAL DEMAND ; 

END; 


PUT EDIT. (SUMGOST (X47 sige a). 
TOTAL AVG MSRT = TOTAL MSRT / SUM DEMAND; 
SUMCOST2 (LAMBDA) =SUMCOST; 


BOD2=0.0; 
LAMBDA = LAMBDA +1; 


IF SUMCOST <BUDGET 
THEN DO; 
OIE (0) <r 
OK =0.5"(OKl -0ke 
END: 
ELSE DO; 
QK2 = Ok. 
OK =0.5*(OK1"> Oka 
END; 
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END; 
Pier SDE T( MOT AVERAGE MSRT IS’, 
bOtdleaveemsceneees ),x(2),F(10,5)); 


END BBO; “DATA 
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APPENDIX B 
COMPUTER PROGRAM FOR N.P.S. MODEL BY PL/I 


*PL/C ATR SOURCE L=4000 P=70 T=(3,30) XREF 
OPTIMAL : PROCEDURE OPTIONS (MAIN) ; 
DCL (ITEM, ITEMS, BASE, BASES; 1, X, Ky BASE@ot eq cee 
depot _ stock,OPTI_ DEPOT STOCK(6),SUM,TOTAL STOCK) Fixaie 
DCL (COST(9,9,6),BASE REP TIME(9,9,6),BASE DEMAND(9,9,6), 
~ ORDER _TIME(9,15,6),dEPOT REP TIME(9,15,6), 
READY RATE(5),TIJ,NEW PROB,base prob(9,15,6))FLOAT; 
DCL (DEPOT DEMAND,SUM DEMAND ,MU,Y3,PX, POISSON PROEe Une 
msrt ,SUM9 ,SUM10,,BASE MU,SUM BACKORDER, BASE BACKORDER, 
TOTAL COST,LAMBDA T,EXPECTED BACKORDER) FLOAT ; 
DCL (TOTAL DEMAND, TOTAL MSRT,TOTAL AVG MSRT,item cost, 
base res _time(/))float; 
LEEMS == —s- 
BASES = 3; 
TOTAL COST =0.0; 
DO Lf =1°70O TEMS, 
DO J il TOSS aS ns: 
Koes es 
GET LIST(COST(1,J,K),BASE REP SEIME ti oe 
BASE DEMAND(1,J,K+2) ,ORDER TIME (hr 
BASE PROB(1,J,K+4) ,DEPOT REP TIME(1 (3K -ooem 


END; 
Ee 
TOTAL DEMAND = 0.0; 
TOTART Moki = OU 
DO: LREMe == tO ee MSE 

SUM = 30 

DEPOT DEMAND = 

SUM DEMAND = QO. 


O02 
0 


>i 
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DO BASE =1 TO BASES; 
TIJ =(1.0-BASE PROB(ITEM,BASE,5))*90 + 
BASE PROB(ITEM,BASE,5)*BASE REP_TIME(ITEM, BASE, 2); 
MU =TIJ*BASE_DEMAND(ITEM, BASE, 3); 
RUecKiP intone Megers, ',TIJ): 
PUT SKIP LIST('BASE MU WHEN TIJ DAYS',MU); 
NEW PROB = .74; | 


PUT SKIP LIST('NEW PROBILITY IS‘ ,NEW PROB); 
Keel; : | . 
Wig 1. 0S 
POISSON = EXP(-MU); 
PROB SUM = POISSON; 
DO WHILE (PROB SUM <= NEW PROB); 
Neo, = eo eK : 
PROBE == TOlSSON = Y3: 
PROB SUM = PROB SUM + PROB; 
Ke ee 
END: 
BASE STOCK(ITEM,BASE) = K -l; 


SUM =e othe eBASE STOCK (LITEM, BASE); 
DEPOT DEMAND =DEPOT DEMAND + BASE DEMAND 
(item,base,3)*(1.0 -BASE PROB(ITEM,BASE,5)); 
SUM_DEMAND = SUM DEMAND + BASE DEMAND(ITEM,BASE,3); 
END; 


DEPOT STOCK s=0;. 
Mmern = 5.5; 
DO WHILE (MSRT >= 5.2); 
LAMBDA _T = DEPOT DEMAND * DEPOT REP TIME 
(item,bases,6); 
PX = EXP(-LAMBDA T); 
SUM10 = LAMBDA T - DEPOT STOCK; 
IF DEPOT STOCK > 0 


bY 


THEN =DO'; 
DO X = 0 TO DEPOT SStOce 
SUM10 = SUMIO + (DEPOT STOCK N= a. 
PX = LAMBDA 1 ~]ie) eer 
END; 
END; 
EXPECTED BACKORDER = SUM1Q; 
DEPOT DELAY =EXPECTED BACKORDER / DEPOT DEMAND; 
SUM BACKORDER = 0.0; 


DO BASE = 1 TO BASES; 

READY RATE(BASE)=0.0; 

BASE RES TIME(BASE) =BASE REP TIME(ITEM, BASE, 2)* 

. .,BASE PROB(ITEM,BASE,5) + (l- | 

BASE PROB(ITEM,BASE,5))*(ORDER TIME(ITEM, 
BASE,4)+DEPOT DELAY) ; 

BASE MU = BASE RES TIME(BASE) * BASE DEMAND 

(item,base,3); 

PUT SKIP LIST('BASE MU IS',BASE MU); 

PX = EXP(-BASE MU); 

SUM9 =BASE MU - BASE STOCK(ITEM, BASE); 


IF BASE STOCK (ITEM; BASE) 0 
THEN TDG; 
DO X =0 TO BASE STOCK( ITEM, BASE) ; 
SUMS = SUM9 +(BASE STOCK(ITEM, BASE) ame 
PX; 
READY RATE(BASE) = READY RATE(BASE)+PX; 
PX = BASE MU ~ PXO Ouse 
END; 
END; 
BASE BACKORDER = SUM9; 
SUM BACKORDER = SUM BACKORDER + BASE BACKORDER; 
END; 
MSRT =SUM_ BACKORDER /SUM_ DEMAND; 
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PUT SKIP LIST('MSRT IS‘ ,MSRT); 

DEPOT STOCK = DEPOT STOCK +1; 
END; 
TOTAL DEMAND = TOTAL DEMAND + SUM_DEMAND; 
TOTAL MSRT = TOTAL MSRT + MSRT*SUM_DEMAND; 


Orie OL ws tOCK(fthM) — DEPOT STOCK -1; 
oP eolOCkwe——scUMerOrTl DEPOT STOCK(ITEM) ; 
iiemecdsl = [Or slOck.~ COST(ITEM,BASES, 1); 


Mileowir PoE Ohll DEPOT STOCK(ITEM) )(X(2),F(6)); 
DO BASE = 1 TO BASES; 
. PUT EDIT(BASE STOCK(ITEM,BASE) ,READY RATE(BASE) ) 
Ca Zee) 5 G26. 4) ); 
END; 
Bees it (aor) (X(2),F(9,4)); 
PULSOKIP (3 )5 
pOrAL COST s=—TOTAL COST + ITEM COST; 
END; 
Pio PIT, TOTAL COST IS ,TOTAL COST)(X(3),A(14), 
22) 
TOTAL AVG_MSRT = TOTAL MSRT / TOTAL DEMAND; 
Pobeokib EbIl( TOTAL AVERAGE MSRT IS’ ,TOTAL AVG MSRT) 
eee (20), X(2),F(955)); 
END OPTIMAL; *DATA 
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